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Abstract 



The /3-function is calculated for an SU(N) Yang-Mills theory from an ansatz for the 
vacuum wavefunctional. Direct comparison is made with the results of calculations of the 
/3-function of QCD. In both cases the theories are asymptotically free. The only difference 
being in the numerical coefficient of the /3-function, which is found to be —4 from the 
ansatz and —4+ ~ from other QCD calculations. This is because, due to the constraint of 
Gauss' law applied to the wavefunctional, transverse gluons (which contribute the |) are 

d ' 

omitted. The renormalisation procedure is understood in terms of 'tadpole' and 'horse- 
shoe' Feynman diagrams which must be interpreted with a non-local propagator. 



* On leave of absence from ITEP, B.Cheremyshkinskaya 25, Moscow, 117259, Russia. 
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1. Introduction. 



One of the main problems in modern quantum field theory is the understanding of 
low energy phenomena in QCD, such as confinement, chiral symmetry breaking, or, in 
more general terms, the strong coupling problem. To have analytic results for the ground 
state of an asymptotically free non-abelian gauge theory, with the associated enhanced 
understanding of the underlying physics, would be invaluable in the understanding of 
these phenomena. Although many promising ideas have been suggested in the first nearly 
quarter century of QCD, e.g. 0], we are still far from a completely satisfactory answer. 

The arsenal of non-perturbative methods available to tackle strongly interacting con- 
tinuum theories is limited. Although in simple quantum mechanical problems a variational 
approach is often easy to use - it is usually enough to know a few simple qualitative fea- 
tures in order to set up a variational ansatz that will give good results for the ground 
state energy and other vacuum expectation values - this is not the case in quantum field 
theories, the complexities of which pose difficult problems, as discussed by Feynman, || . 

Recently, a variational approach based on a gauge invariant Gaussian wavefunctional 
has been studied. This method was applied to QCD and QED 3 , H and j4| respectively, 
where, although in its infancy, it has independently verified many old, and given some 
new, results. A Gaussian approach to the wavefunctional of QCD (the so-called squeezed 
gluons) has been studied in many different papers, see || and the references therein. In 
the case of QED3 the variational method reproduced Polyakov's path integral results of 
the mass gap and string tension of the theory, 

The variational calculation carried out for the SU(N) purely Yang- Mills theory in 3 + 1 
dimensions, ||, has found that the ground state energy is minimal for a state which is 
different from the perturbative vacuum even though the perturbative vacuum state was 
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included in the variational ansatz. Dynamical scale generation takes place and the gluon 
(SVZ) condensate in the best variational state was found to be non-zero. 

The ansatz used for the vacuum wavefunctional was Gaussian but it was also required 
to be gauge invariant. To satisfy this, the Gaussian wavefunctional is projected onto the 
gauge invariant sector. The form of the wavefunctional is discussed briefly in the next 
section and in detail in || . The variational parameter of the theory is the mass scale, M. 
A non-zero value of M for the minimal state corresponds to a non-perturbative dynamical 
scale generation. 

Within the variational calculation of || it was conjectured that the variational ansatz 
proposed yields a theory with a coupling constant that runs as the coupling constant of 
(asymptotically free) QCD. In this paper we prove this conjecture. It is shown that the 
proposed variational ansatz yields an effective non-local, non-linear sigma model in three 
dimensions which, when renormalised to first order, has practically (a precise qualification 
of this word will be given below) the same /^-function as asymptotically free QCD. Further, 
the renormalisation procedure is interpreted in terms of the Feynman diagrams included 
and it is found that, due to the non-local nature of the propagator, a new 'horse-shoe' 
diagram is non-zero and makes a vital contribution. 



In this paper, the /3-function for the charge of the variational ansatz is calculated to 

be 



rif is the number of species of fermions in representation r, which is zero for the purely 
Yang-Mills model. For SU(N), 02(G), the quadratic Casimir in the adjoint representation, 
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is N. We immediately see that the only difference between the two /3-functions is in the 
numerical terms 4 and (4 — ~). The 4 is due to the ant i- screening effect of longitudinal 
gluons. The | is due to the screening effect of virtual transverse gluons. Identically 
the same /3-function as QCD is not obtained because only longitudinal gluons have been 
included within the ansatz. The longitudinal gluons are included because they satisfy 
Gauss' law with a source, which is used as a constraint upon the vacuum wavefunctional 
in setting up the variational ansatz. 

In this paper the renormalisation calculation to first order is presented and the renor- 
malisation of the effective charge in the variational wavefunctional is obtained. In the 
first section we shall describe the variational ansatz of |§ and show how this leads to 
an effective non-local, non-linear sigma model in three dimensions. In the second section 
we perform the renormalisation group transformation, integrating over high momentum 
dependent modes to yield an effective action for the low momentum modes with a renor- 
malised coupling constant, and we interpret the renormalisation procedure by considering 
the Feynman diagrams which contribute. 



2. The variational ansatz. 



For a full discussion of the variational ansatz and all details of the subsequent vari- 
ational calculation the reader is directed to the original paper, ||. An overview of the 
variational ansatz and the form of the effective action is given in this section. 



The SU(N) gauge theory is described by the Hamiltonian, 

H = I d 3 x \ l -Ef + l -Bf 



(2.1) 



where 
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E?(x) = i 



5At(x) 
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B?(x) = -e ljk {d,At{x)-d k A«{x)+gr bc A%x)Al{x)} (2.2) 
and all physical states must satisfy the constraint of gauge invariance (Gauss' law); 

G a (x)V[A] = [diE?(x) - gf^A^E^x)] V[A] = (2.3) 
Under a gauge transformation U (generated by G a (x)) the vector potential transforms as 

A«(x) -> Af a (x) = S a \x)A\{x) + A?(x) (2.4) 

where 

S ab ( x ) = -tr (r a WT b u) ; X^x) = -tr (r^t^f/) (2.5) 

and r a are traceless Hermitian N by N matrices satisfying tr(r a r b ) = 25 ah . For SU(3) the 
algebra and structure constants are defined, for example, in ||. 

The initial ansatz for the ground state wavefunctional is of vital importance. It must 
incorporate the properties of all such physical states and yet it must not lead to a solution 
which is incalculable if any progress is to be made. In this formalism one calculates 
expectation values of local operators with the ansatz for the ground state, 

< O >= - [ D<f)^*[(f))0^[(f>] (2.6) 

and then minimises with respect to the variational parameter. A calculation of this 
kind is tantamount to evaluation of a Euclidean path integral with the square of the 
wavefunctional playing the role of the partition function. One should therefore be able to 
solve exactly a <i-dimensional field theory with the action 

S[ct>} = -log¥*M*M (2.7) 



Since in dimension d > 1 the only theories one can solve exactly are free field theories, 
the requirement of calculability almost unavoidably restricts the possible form of the 
wavefunctional to a Gaussian (or, as it is sometimes called, squeezed) state: 

*[A*} = exp {-1 J d*xd 3 y [A?(x) - Q(x)] (G^(x,y) [A*(y) - <*(y)]} (2.8) 

with ((x) and G(x,y) being c-number functions. The requirement of translational invari- 
ance usually gives further restrictions: ((x) = const, G(x,y) = G(x — y). 

There is, however, one obvious difficulty with this idea. It is impossible to write down 
a Gaussian wavefunctional which satisfies the constraint of gauge invariance. Under the 
gauge transformation the wavefunctional transforms as 

*[A a ] - 9[{A u )t] (2.9) 

In the abelian case it is enough to take diG^ 1 = to satisfy the constraint of gauge 
invariance. In the non-abelian case, however, due to the homogeneous piece in the gauge 
transformation (|2.4j) , no gauge invariant Gaussian wavefunctional exists. 

The proposed solution to this problem was to simply project the Gaussian wavefunc- 
tional onto the gauge invariant sector and to restrict the calculation to the case of zero 
classical fields (( = 0). The variational ansatz is therefore 

V[A1] = J DU(x) exp |-1 J d 3 xd 3 y Af a (x)G^(x - y) Af{y)} (2.10) 

with A^ a defined in ( |2.4| ) and the integration performed over the space of special unitary 
matrices with the SU(N) group invariant measure. 

Further restrictions upon the form of G lead to considerable simplifications. Firstly, 
only matrices of the form 

G? j {x-y) = 8 ah 8 ij G{x-y) (2.11) 
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are considered. This is certainly the correct form in the perturbative regime. If it was not 
for the integration over the group, G^ 1 would be precisely the (equal time) propagator of 
the electric field. Due to the integration over the group, however, the actual propagator 
is the transverse part of G^ 1 . The longitudinal part diGjl drops out of all physical 
quantities, giving, without any loss of generality at the perturbative level, Gij ~ 5y. 
Also, in the leading order in perturbation theory, the non-abelian character of the gauge 
group is not important. The 5 ab structure is then obvious. 

The form of G can be restricted further using additional perturbative information. 
The theory is asymptotically free. This means that the short distance asymptotics of 
correlation functions must be the same as in the perturbation theory. Since G' 1 is directly 
related to correlation functions of gauge invariant quantities in perturbation theory, it is 
taken to have the form, 

G-\x)^- 4) x-^0 (2.12) 

The non-perturbative theory is also expected to have a gap. In other words, the correlation 
functions should decay to zero at some distance scale, 

G{x)~0, x>^ (2.13) 

The variational ansatz is built in the simplest possible way. M is taken to be the only 
variational parameter and this is done by choosing G(x) to be of a form that has the ultra- 
violet and infra-red asymptotics of (|2.12|) and ( [2.13|) . A non-zero result for M means a 
non-perturbative dynamical scale generation in the Yang-Mills vacuum. The form of G^ 1 
used in this paper and the variational calculation of || has the Fourier transform 

rt)J^ iffc 2 >M 2 

U [ > 1 M iik 2 <M 2 1 > 



Equation (|2.10 ) together with equations (|2.11 ) and ( |2.14 ) define our variational ansatz. 
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It will now be shown that the action (|2.7|) is in fact a non-local, non-linear sigma model. 
Again, the reader is referred to || for the original account. With the given ansatz, the 
remaining problem is to calculate expectation values of local operators, such as, 

< O >= - / DUDU' <0> A (2.15) 

Zj J 



where 



< O > A = J DAe~* J ' dxd V AY a ( X )G- l ^-y)Af a (y) q e -\ J dx'dy'Af b (x')G- 1 (.x'~y')Af b (y') 

DAe~^ I dxd y A Y a ( x ) G ~ 1 ^-y)AY a (y)Oe'^ I dx , d y 'A b j (x')G- 1 (x'- y ')A b j (y') ^.16) 



where, since only gauge invariant operators are to be considered, the change of variable 
Ai — > AY U has rendered one of the group integrations redundant. 

For convenience, the definition 

a?(x)= I <Pyd z z\\{y)G- l {y - z)S bc {z)(M- 1 ) ca {z,x) (2.17) 



is made so that the Gaussian integration over Ai is / DAexp[—^(A+a)A4(A+a)]. Chang- 
ing variables and performing the integration yields the following form of the normalization 
factor Z, 

Z = J DU exp{-T[U}} (2.18) 

with an effective action 

F[U] = -TrlnX + -A a A ac A c (2.19) 

where 

A ac {x, y) = [G{x - y)5 ac + S a \x)G(x - y)S Tbc {y)]- 1 (2.20) 

is the 'effective inverse propagator' and multiplication is understood as the matrix multi- 
plication with indices: colour a, space % and position (the values of space coordinates) x, 
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i.e. 

(AB)%(x,z) = J d 3 yA« b (x,y)B bc k (y } z) } XOX = J d 3 xd 3 y\^x)0%(x - y)X)(y) 

(2.21) 

The trace Tr is understood as a trace over all three types of indices. In equation ( |2.19| ) 
we have defined 

S$(x,y) = S ab (x)5 %] 5(x - y), M$(x,y) = [S Tac (x)S cb (y) + 5 ab ]G-\x - y)5 ij (2.22) 

where S ab (x) = |tr (r a UV b u) and A? (x) = |tr (r a WdiU^ were defined in fl2j) and tr 
is a trace over colour indices only. One should also note here another useful definition, 
the completeness condition for SU(N); 

t^i = 2(M* - ( 2 - 23 ) 

The path integral ( |2.18| ) defines a partition function of a non-linear sigma model in 
three dimensional Euclidean space. The action of this sigma model is rather complicated. 
It is a non-local and a non-polynomial functional of U(x). We shall see how the coupling 
appears in the effective action in both high and low momentum cases. Various approxi- 
mations are made here which are explored much more rigorously in the calculation of the 
succeeding section. 

For high momentum modes, with the standard parametrization U(x) = exp[i|0 a (x)r a ], 
one gets A°(x) = — di<f) a (x) + 0(g), S ab (x) = 5 ab + 0(g) and the leading order term in the 
action becomes: 

^ J d 3 xd 3 yd t r(x)G-\x - y)d l( t> a (y) (2.24) 
This is just a free theory with a non-standard propagator; 

<nx)<p b (y)>=25 ab [d!d!G-\x-y)r=25 ab J ( g 3 eXp[ ^ = y)] (2.25) 
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For the low momentum modes, to a first level of approximation, the space dependence 
of S (x) is ignored in the term SGS T giving, with the fact that S is an orthogonal matrix, 
the approximation SGS T — > G. Then, using the completeness condition, ( j2.23|) , and the 
fact that tr(Wd t U) = we can write 

A?(x)A?(x) = -(l/# 2 )tr (ftfdiU) tr (r'rfdiU) = -{2/g 2 )tx (Wd l U U^diU) (2.26) 

In this approximation the action becomes 

1 t Mr 
— tr J d 3 xd 3 yd i U^x)G-\x - y)diU(y) = — tr J d 3 xdiU^x)diU(x) + ... (2.27) 

where +... corresponds to all the higher terms in g. 



3. Calculation of the /^-function. 



Having obtained the effective action of ( |2.19 ), which is a non-local, non-linear sigma 



model in three dimensions, we shall proceed to calculate the j3 function. The form of the 
/3-function is deduced from the renormalised coupling constant, in the spirit of |7| where 
a similar calculation was performed for the non-linear sigma model in two dimensions. 
We shall perform the renormalisation group transformation by integrating over high mo- 
mentum dependent modes (containing Fourier components k > M) leaving an effective 
action for low momentum dependent modes (containing Fourier components k < M). 
The coupling constant of the effective action is renormalised to first order; up to terms 
quadratic in the high momentum modes. A physical interpretation of these terms by con- 
sidering the corresponding Feynman diagrams is given in section 4. The high momentum 
modes are the quantum field, and the low momentum modes are the classical field, of the 
background field method. 

10 



Next we shall discuss the decomposition of the group elements into high and low 
momentum dependent modes and in the second subsection we shall explicitly calculate 
the j3-f unction. 

3.1 Quadratic approximation for the high momentum modes. 

The ansatz proposed for the decomposition of group elements into high and low momen- 
tum dependent modes is 

U(x) = U L (x)U H (x) (3.1) 

where Ul(x) contains Fourier components k < M and Uh{x) contains Fourier components 
k > M. This can be considered as a decomposition of the group parameter, <fi a (x). 
If we write U(x) = exp[^<f) a (x)T a ] then the decomposition can be written as <ft a (x) = 
4>h( x ) + 4>l(. x )i where (j) a H L (x) are the group parameters of Uh,l{x), respectively, and (j> a H 
and 4> a L are taken to be orthogonal; <p a H (x)(p a L (x) = 0. 

Written in terms of the group parameters, we can explicitly see that this ansatz has 
similarities to that used by Polyakov in his treatment of the non-linear sigma model in 
two dimensions, 0, but without the normalisation of 4>%(x) constrained to be that of 
<j) a (x). In his ansatz, Polyakov forced \4> a \ 2 = \4>1\ 2 i proposing 

= a L (x)(l - |^| 2 )5 + 0^(x)ef(x) (3.2) 

which has been written in the notation of this paper. e"(x) form a complete basis of unit 
vectors orthonormal to 4>1(x). We do not use this ansatz for the decomposition, however, 
because for each order of g in the calculation it mixes the high and low momentum group 
parameters making the desired decomposition (|3.1|) difficult. 

For the rest of this section we shall employ the decomposition ( |3.1| ) to write the 
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effective Lagrangian for A" L up to terms of 0(g 2 ). This corresponds to all terms up to, 
and including, those quadratic in the field 0. 

It is necessary to first note two identities for A" and S ab . Using the completenes 
condition for SU(N), (|2.23|), with the cyclic property of the trace and the fact that 



tr[Ur a U+] = tr[r a ] = we find, 

S ac ( x ) = S%(x)S b L c (x) (3.3) 

and 

= S^Xl L {x) + Xl H {x) (3.4) 

Using the same mathematical properties, we should also note that S ab is an orthogonal 
matrix, 

S ab (x)S Tbc (x) = S ab (x)S cb (x) = ^tUU^U^x^U^U^x))^ = \tr[r a r c } = 5 ac 

(3.5) 

First, we shall evaluate the inverse effective propagator, A ac (x, y) (|2.20 ), up to terms 
quadratic in the coupling constant. To do this we write, 

U H (x) = eX p[|^(x)r a ] = 1 + % ^\x)r a - 9 ^{<f{x)r a f + 0(g 3 ) (3.6) 

We shall use the following normalisation of the generators of SU(N), (e.g. with constants 
defined for SU(3) in §), 

[r a ,T b ] = 2if abc r c (3.7) 
Itr[rV 6 ] = 5 ab (3.8) 

We find, 

S a H d (x) = \tr[T\T d ]-^<p\x)tr[T a {T b y\] (3.9) 
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-^ ( p b ( x )^( x )tr[T a [r b , [r c , r d }}} + 0(g 3 ) 
lb 

= S ad - gf adb <p\x) - £f<*f**<f(x)<lf(x) + 0(g 3 ) 

Therefore, using (|3.3|) and considering the low momentum group elements, Ul(x), to be 
slowly varying, such that UjJ^U^y) — 1 and S'^'(x)Sf } (y) ~ 5 de , we find, 

S ab (x)S Tbc (y) ~ S a H d (x)S cd (y) (3.10) 
= 5 a *-gA™(x,y)-^A™(x,y) + 0(g 3 ) 

where 

AT(x,y) = r b {<?{x) - <j> h {y)) (3.11) 

A?(x, y) = f^f^ix^ix) + <p{y)P{y) - 2<f> b (x)^(y)) (3.12) 



We now re-write A ac (x,y) as, 



A ac (x,y) = G-\x-y)[28™\-\l-g [ —^Al c {x,y) (3.13) 

9 2 i 26a r\ A a C ^ y) + 0{g 3 rl 



2 tr[5 a 

sjac 2 
= Y G ( X - V)[1 + ^^} Ar ^ V)]+0{93) 

where [c^]" 1 = 5 ac , 5 ac [5 ac ]~ 1 = tr[5 aa ]. Also, 

A?(x,y) = (3.14) 
A?(x,y) = C 2 (G)(0 fe (a;)0 6 (a;)+0 6 ( 2/ )0 6 ( 2/ )-20 b (a;)0 6 ( 2 /)) (3.15) 

C2{G) is the second Casimir operator and G denotes the adjoint representation in this 



case. For SU(N), C 2 (G) = N. Equation ( |3.15|) is obtained using the relation, 

jacdfbcd = c 2 (G)5°* (3.16) 
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Next we need to evaluate Xf(x) up to 0(g 2 ). Using (|3.4j) , (3.9) and following similar 
intermediary steps to the calculation of (3.9), we find, 

KM X ) = -di<f> a (x)+0(g 3 ) (3.17) 

A?(x) = \l L {x)-d l ct ) \x)-gr b % L {xW^) (3-18) 
- y f ace f bde ^i,L ( x ) 4> c { x ) d ( x ) + 0(g 3 ) 

Therefore, we can now write the effective action as, 

T[U} = J d 3 xd 3 yT[x jy ] (3.19) 



where 



V{x,y) = - 4 d4 a {x)G-\x-y)d i r{y) (3.20) 
+hl L ( x )G-\x-y)Xl L (y) 



+ g -F b9 i\ b l AxW(x)G-\x - y )d4 a {y) 

+d i r(x)G- 1 (x-y)\l L (y)r(y)} 

_£ f a 9 e fbd e xbiL{x)G -l {x _ y)^ L (y)[^(x)^(x) 

+<f> 9 (y)Ay) - 2^(x)/(y)] 
+ 9 y Q \Ux)G-\x - Z/)A^(y)^|[0 b (x)0 fe (x) 
- 2<p b (x)<f ) b (y)] 

where T(x,y) is the effective Lagrangian. 

3.2 Renormalisation group transformation. 

The renormalisation group transformation is performed by integrating over the high mo- 
mentum dependent field, <fi a (x). This is akin to integrating over the fluctuating (quantum) 
fields in the background field method. 
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We consider 



T°(x,y) = - A d4 a {x)G- 1 {x-y)d i 4> a {y) (3.21) 



to be the zeroth order Lagrangian for the field 4> a (x). All other terms involving <fr a are 
treated as perturbations of this Lagrangian. ( ^.21[ ) implies the two-point Green's function 

< n^iv) >= 25 ab {d*d!G-\x - y)]- 1 (3.22) 



G 1 (x — y) is defined in ( |2.14 ). So, because (p a (x) is defined to have Fourier components 
k > M, we can write, 

<^)/ W>=2 ^/g^ gPlto)] ,3.23) 

with the integral performed over the limits M < k < A, < <p < 2tt, < 9 < it, where A 
is the ultra-violet cut-off. 

Treating terms other than T as perturbations, we see that 

J D<f>esxp[ J d 3 xd 3 y{T°{x,y) + = exp[ J d 3 xd 3 y < F[<f>] >} (3.24) 

So we now write the effective Lagrangian as 

rz(*,y) = \\t L (x)G- 1 (x-y)\l L (y) (3.25) 
+ \f ah9 [\l L (x)G~\x -y)< (j> 9 {x)d i( j> a {y) > 
+ < drf* (x)^(y) > G-\x - y)Xl L {y)} 
-j r e f bde \l L {x)G-\x-y)\l L {y){< <P 9 {x)ct> d {x) > 
+ < ^{y)ct> d {y) > -2 < 9 (x)0 d (y) >] 
+ ^KL^)G-\x-y)Xl L (y)^l^[< <j>\x)<p\x) > 
+ < <j> b (y)<l> b (y) > -2 < <j> b {x)<j> h {y) >] 
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x y 
Fig. 1 

First considering the terms of 0(g) we see that both their contributions are zero. This 
is explicitly seen as both of the correlations < <p 9 (x)di<j) a (y) > and < di<p a (x)<f) 9 (y) > have 
the structure 5 a9 preceeded by the totally antisymmetric group structure constant f abg . 

Considering the terms of 0(g 2 ), we need to examine the evaluation of ( |3.23| ) in some 
detail. First, we note that, 



fi ab A 

< >=< <t> a (y)4> b (y) >= — log — (3.26) 

7T Z M 



which occurs in ( |3.26|) in terms such as 



9 ^XUx)G-\x - y )A^(y)^M < 4> b (x)4> b (x) > (3.27) 

This can be represented by a Feynman diagram, Fig 1, which shows a tadpole dia- 
gram. The external lines correspond to low momentum fields, U(x), and the internal 
loop represents the integration over high momentum fields. The dotted line corresponds 
to G^ 1 (x — y). It is important to note that in the region corresponding to k < M, 
G^ 1 (x — y) = M5(x — y) and the propagator becomes local, associating either end of the 
dotted line, and the standard tadpole diagram is recovered. The evaluation of fl3.26p is, 
however, the same at all scales of spatial separation, \x — y\. 

To evaluate the other terms in Yi,(x,y), we write (|3.23 ) in the form, 



Kab F A\x—y\ s j n + 

<rM>= -J M ^ dt — (3.28) 
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where the change of variable t — k\x — y\ has been made, t is a dimensionless variable. 
We shall introduce a scale, /i, such that for t < fi, sint ~ t. This allows us to write, 

< <p a {x)<p\y) >= { - . ° n ^\ X ~ V \ >fi (3.29) 

W V ; 1 -^rbg[M|a:-j/|] + 6 M\x - y\ < \i 1 ; 

where a and b are finite contributions which are independent of M and A and hence are 
ignored in the following. All /i dependence is in the finite terms. 

The correlations < <fi a (x)(f) h (y) > appear in Tl(x, y) in terms such as, 

- 9 jXUx)G-\x - y )Ak(y)£MS < ^(x^iy) > (3.30) 

which can also be interpreted in terms of Feynman diagrams. Fig. 2 shows the diagram 
corresponding to (|3.30| ). It depicts a horse-shoe. As in Fig. 1, the external lines correspond 



to the low momentum fields, Ul(x), the horse-shoe represents the integration over the high 
momentum fields and the dotted line again represents G^ 1 {x—y). The horse-shoe diagram 
is unimportant in the region j| < \x — y\ as it only gives a finite contribution, independent 
of a cut-off or M. In the region < \x — y\ < -fa the horse-shoe diagram is important 



and contributes ^ log M ^_ y ^ , ( |3.29| ) . We must not here immediately interpret \x — y\ to 



therefore be the inverse of the cut-off but rather take our interpretation from the form of 
the resulting Lagrangian. 

Substitution of ( |3.26| ) and (|3.29| ) into Tl(x, y) yields the effective action for the low 



momentum modes with a renormalised coupling constant, g. With \°[ L = ^ti[T a U^ diU l\ 
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this can now be written as, 



T L (x,y) = \\l L (x)G- 1 (x-y)\l L (y) (3.31) 



-Lti[d l U L (x)G~ 1 (x-y)d l Ul(y)} 



where 



~ a 2 = \ g 2 (I + ^C 2 (G) log ±) M\x-y\>p 

\g 2 {l + ^C 2 {G)\ogk\x-y\) M\x - y\ < // l ' ' 

(3{g) is calculated using the standard definition, 

m = M m §la ' A (3 - 33) 

but first we need to interpret the appearance of \x — y\ in ( |3.32| ). For \x — y\ > M 



is the renormalisation scale, whereas, for \x — y\ < jj, r^z^\ should be interpreted as the 
renormalisation scale. Therefore, we can re-write g as, 

~9 = 9+ (£^(G0 lo S ^7 + °(9 5 ) (3.34) 



where 



, ,, [ M \x-y\ > & . 

M = < _j_ _ y \ f_ (3.35) 

{ \x-y\ \ X V\ < M 



Thus we obtain, 



r\ 3 

(3(g) = M'—~g\ gA = - * 4<7 2 (G) + 0(g 5 ) (3.36) 
This should be compared to the standard /5-function for QCD, e.g. @, 

M = -(J^K 4 - \)<b{G) - §"/] + 0(g 5 ) (3.37) 

In this paper and in ||, a variational approach to QCD has been considered in the absence 
of fermions, n/ = 0. With a gauge group SU(N), C^G) = iV. The discrepancy between 4 
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in the result of this paper and the y of the standard QCD /3-function quoted above is due 
to the fact that only longitudinal gluons are included in the variational ansatz and that 
transverse gluons, which would contribute the extra |, are omitted. This is as expected 
because in the creation of the ansatz it was required that the wavefunctional obey Gauss' 
law, which is indeed satisfied by longitudinal photons with a source. 

The representation of a /3-function coefficient as the sum of two contributions propor- 



tional to —4 and 1/3 is not new at all. It has been known for a long time, [IU| and ['I I 
that in the background field method the contribution of charged particles with spin S to 
the one-loop /^-function is given by 

(-1) 2S 1 
M) = VU-[(2S) J - 3) (3-38) 

where we omit the group factor. The asymptotically free (for integer spins) "spin" factor 
4S 2 , where 5=1 for the vector field A, gives 4 which is precisely what we have obtained in 
our calculations. As was discussed in ]TT[ , this spin factor is related to the influence of the 
background field on the electric dipole moment density, which includes contributions from 
the time-like polarization states of the massless vector fields in the case of the Feynman 
gauge. In a Coulomb gauge, where the time-like modes do not exist as dynamical degrees 
of freedom, one can see that the electric dipole density appears because of Gauss' law. 
This explains why in our method, where the Gauss' law is implemented by construction, 
we obtained the same result. 

It is interesting to note that the same decomposition 4 — 1/3 comes from the results 
for the calculation of the pre-exponential factor, or the renormalisation of the charge, for 
the BPST instanton within the path integral formalism of QCD, ||12|| . In this work, the 
vector field is split into components, = A a ^ inst ^ + a^, and the action expanded in terms 
of the deviation a° from the instanton field A a ( mst \ Analysis of the resulting path integral 
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and examination of the contribution of the zero-frequency modes yields, 
< 0|0t >fns f d 4 xdp 



const, f -JL£.S$exp[-S + 8 log Mp + $ x ] (3.39) 



< 0|0 T > p . th . J p 5 

where p. th. refers to perturbation theory, |0 >t is the vacuum after time T, d A x is the 
measure of integration over the four coordinates of the centre of the instanton, p is the 
scale of the instanton, M is the introduced cut-off parameter and So = ^r- 3>i denotes 
the contribution of the positive frequency modes. In the limit Mp » 1, $x was evaluated 
to the one loop level by means of ordinary perturbation theory to be 

$! = | log Mp (3.40) 

Substituting (|3.40|) into the argument of the exponential in ( p.39|) , the result of the renor- 
malised charge is obtained, 

^f, = ^-2(A- 1 -)\ogMp (3.41) 

The 4 explicitly came from the evaluation of the zero frequency modes and the — | from 
the evaluation of the transverse positive frequency modes. 



4. Conclusion 

The variational ansatz of has been studied and the renormalisation of the effec- 
tive charge in the wavefunctional has been carried out up to 0(g 2 ). In this procedure, 
the effective action is a non-local, non-linear sigma model in three dimensions where the 
fields considered are the group elements of the original gauge transformation. The group 
elements are decomposed into low and high momentum dependent components and the 
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renormalisation transformation is effected by integrating over the high momentum depen- 
dent modes up to (quadratic) terms of 0(g 2 ). The (3 function is found to be 

P(g) = -j£^^C 2 (G)+0(g 5 ) (4.1) 

This should be directly compared with the known (3 function for QCD, ||; 

Piff) = -(|^[( 4 " \)°^ G ) ~ l n f C (r)\ + 0(g 5 ) (4.2) 

The only difference between the two (in the absence of fermions n/ = 0) is the inclusion of 
the factor | in the latter, which is due to the screening effect of virtual tranverse gluons. 
Only longitudinal gluons were incorporated in the ansatz as they satisfy Gauss' law, the 
constraint ensuring gauge invariance of the wavefunctional. 

The interpretation of the renormalisation procedure by considering the Feynman dia- 
grams shows that a vital contribution is made by tadpole diagrams and new 'horse-shoe' 
diagrams which must be interpreted with the non-local propagator G _1 (x — y). The 
renormalistion scale is found to be M when \x — y\ > 4i and r^— r when \x — y\ < 4r, 
where M is a dynamically generated mass scale (and we have taken /i, a dimensionless 
constant from the previous section, to be 1 for simplicity). The explicit calculations made 
in this paper give a proof of the conjecture made in || that to calculate the gap in a 
variational approach one can use an effective QCD coupling constant and that this gap 
will automatically be related to Aqcd- 
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